Stochastic Reduction in Nonlinear Quantum Mechanics 
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Stochastic extensions of the Schrodinger equation have attracted attention recently as plausible 
models for state reduction in quantum mechanics. Here we formulate a general approach to stochas- 
tic Schrodinger dynamics in the case of a nonlinear state space of the type proposed by Kibble. 
We derive a number of new identities for observables in the nonlinear theory, and establish general 
criteria on the curvature of the state space sufficient to ensure collapse of the wave function. 
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A generalisation of quantum mechanics was considered 
by Mielnik jl], who introduced the notion of nonlinear 
observables. Two alternative extensions of the standard 
quantum theory were then proposed by Kibble ||^. The 
first alternative is based on the phase space formulation 
of quantum mechanics. Here, if we work with the space 
of rays through the origin of Hilbert space, then the 
Schrodinger equation reduces to Hamilton's equation of 
classical mechanics , except that the quantum Hamil- 
tonian is of a special restricted form. Thus a natural 
generalisation is to remove this constraint. When such 
trajectories are lifted from the space of rays to Hilbert 
space, we obtain nonlinear wave equations. 

The general properties of nonlinear observables were 
subsequently analysed in detail by Weinberg Q . Follow- 
ing this, it was pointed out by Gisin |^ that the evolution 
of the density matrix is not autonomous in the nonlin- 
ear mechanics of (|,^, and that this may be physically 
undesirable. However, it was also indicated in that 
there is another type of nonlinear quantum dynamics for 
which the evolution of the density matrix is autonomous, 
and a number of desirable features of linear evolution are 
extended in a natural way. This is the stochastic dynam- 
ics developed by Pearle and others These dynamics 
are of significance because they exhibit natural reductive 
properties: starting from a given initial state, the system 
evolves stochastically in such a way to ensure collapse to 
an eigenstate of one or more designated observables. 

Kibble's second alternative for a nonlinear quantum 
theory is in essence to consider a general Kahler mani- 
fold as the phase space of quantum mechanics, instead 
of the space of rays. The idea is that, in the presence of 
interactions, the states accessible to a quantum system 
constitute a curved space 971 which has the structure of a 
complex manifold endowed with a compatible symplectic 
structure. The dynamics of the state are then governed 
by a Hamiltonian flow which is also an isometry. 

In the present article, we consider stochastic state re- 
duction models within the framework of Kibble's second 
theory. The advantage of a stochastic dynamics in this 
context is that it leads to a probabilistic interpretation, 
a feature hitherto missing in the nonlinear theory. Re- 



markably, many of the key features of the basic stochastic 
reduction models carry through to a fully nonlinear state 
space. The main results are to determine general crite- 
ria sufficient to ensure state reduction in the nonlinear 
theory, and to express these criteria directly in terms of 
geometrical features of the state manifold. Thus it is the 
geometry of the quantum state manifold that determines 
whether reduction takes place, and if so, how rapidly. 

After introducing the relevant state space geometry 
and elements of stochastic calculus on manifolds, a num- 
ber of identities concerning the properties of quantum 
observables are established in Lemmas 1-5. These results 
are then applied to formulate general theorems governing 
reduction processes on nonlinear state spaces. 

Let us first recall briefly the phase space formulation of 
quantum theory. We consider a flnite dimensional com- 
plex Hilbert space TL of which a typical element is de- 
noted (a — 0, 1, • • • , n). Given the Hamiltonian op- 
erator Hp, the dynamics of the state is determined by 
the Schrodinger equation ihdtip" = H^i/;^. The expecta- 
tion Fl^ipaip^ /iJ-ytp'^ of an observable Fg in the state 
is invariant under the scale transformations tp" —> Xip"' 
(A G C — {0}). Hence we can work with the space of 
equivalent classes of state vectors modulo such transfor- 
mations, i.e., the complex projective space P". 

We regard P" as a real manifold F of dimension 2n. It 
is known that F has a natural symplectic structure uJat, as 
well as a Riemannian structure given by the Fubini-Study 
metric gab- These two structures are compatible in the 
sense that there exists an integrable complex structure 



J\ on F, satisfying J^c-^l 



-S% such that VaJ ^ = 



and g°''^ujcb = J%, where Vq is the covariant derivative 
associated with gab, and gacg'^^ = Sj'. We use Roman 
indices (a, 6, • • • ) for tensorial operations on P. The com- 
patibility conditions make P a Kahler manifold. 

The special feature that identifles P as the quantum 
phase space is that the Schrodinger equation can be ex- 
pressed in the Hamiltonian form hdx°'/dt = 2LU°-''VbH{x). 
Here = g°''^g^'^ujcd and H{x) is the Hamiltonian func- 
tion on P, given by the expectation of the operator 
in the equivalence class of state vectors corresponding to 
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the point x G T. The vector = uj°-''\/bH tangent to 
the Schrodinger trajectory satisfies the Kilhng equation 
V(a-Z^f)) = iff H{x) is the expectation of a quantum ob- 
servable in the state x. Therefore, tlie Schrodinger evo- 
lution preserves the distance, and hence the transition 
probability, between any given two states. More gener- 
ally, any isometry of F is a Hamiltonian flow associated 
with a quantum observable. The energy eigenstates are 
the fixed points of the flow, at which Va-ff = 0. The first 
alternative of Kibble is to replace the observable H{x) by 
a general function on F. Then the resulting trajectories 
are Hamiltonian but no longer Killing, and the implied 
dynamics on Ti is governed by a nonlinear wave equation. 
This is not the generalisation we consider here. 

For the consideration of Kibble's second alternative, 
it will be useful first to develop a differential geometric 
framework for the standard operations of quantum me- 
chanics (we set h = 1). U F{x) and G{x) are observables, 
the expectation of their commutator is also an observ- 
able, given by the Poisson bracket 2w°^Vai^Vf,G. Then 
if Xt is a Schrodinger trajectory and Ft = F{xt), it fol- 
lows that dFt = 2uj°'''VaF\/bHdt, where H is the Hamil- 
tonian. This tells us how the expectation of F changes 
along the flow generated by H. For any observable F{x), 
we define the associated dispersion by 

- g'^'VaFVbF (1) 

In the linear theory V^{x) is the squared uncertainty 
of F in the state x. As a consequence of the inequal- 
ity {gabX^X''){gabY^Y'') > (gabX^Y'')^ + 
which holds for all X° and we obtain the Heisen- 
berg relation V^V^ > [w^'^VaFVbGf if we set X"" = 
w"''V&F and Y'' = uj''''VbG. 

Kibble's second alternative for a nonlinear quantum 
theory is to let the state space be a general Kahler man- 
ifold OJt, with metric gab, symplectic structure Wab, and 
complex structure J^. In the nonlinear theory we say 
that a real function F{x) on 371 is an observable iff the 
corresponding Hamiltonian vector field = (jj°'^VbF 
is an isometry. This agrees with the usual characterisa- 
tion of observables when 971 = F. We note that if 971 is 
compact and has vanishing first Bctti number, then any 
Killing field on 971 is Hamiltonian, with at least two dis- 
tinct eigenstates Q . As in the linear theory, we interpret 
F{x) as the expectation of the result of a measurement 
of the given observable in the state x G 971. 

Lemma 1. If F{x) and G{x) are observables, then 
their commutator is an observable. 

The proof is as follows. If X° and Y"" are Hamiltonian 
flows, then their Lie bracket is X^VbF" - r^'VbX" = 
w''''Vbiu}cdX''Y'^). UF{x) andG(a;) are generators of X'' 
and F", then ujcdX^'Y'^ = Lo'"^VcFVdG. Furthermore, if 
X°- and Y°- are Killing, so is their Lie bracket. Therefore, 
the Hamiltonian flow generated by the commutator of 
F{x) and G{x) is Killing. As a consequence we obtain 
also the following nonlinear generalisation of an identity 
due to Adler and Horwitz 



Lemma 2. If F{x) and G{x) are observables, then 
VbFV^VG - VbGV^V'F = uj''''Vbiuj'"^VcFVdG). 

As in the linear theory, the eigenstates of an observable 
F{x) are the points of 971 at which VaF = 0. The value of 
F{x) at a critical point is the corresponding eigenvalue. 
In the nonlinear theory V^{x) does not in general have 
an interpretation as a moment, but nevertheless remains 
a measure of the dispersion of F{x) in the given state. In 
particular, the Heisenberg relation holds. 

The Schrodinger trajectories in the nonlinear theory 
are generated by a Hamiltonian H{x), which we assume 
to be an observable. The following result shows that 
for any observable commuting with the Hamiltonian, its 
dispersion is constant along the Schrodinger trajectory. 
This fact will be used later to derive the stochastic dy- 
namics of the energy dispersion. 

Lemma 3. If F is an observable that commutes with 
the Hamiltonian H, then uj°-^VaHVbV^ = 0. 

The proof is as follows. Equation (jl|) implies that 
uj'''>\/aH\/bV^ = 2w'^^Vai?Vf,VcFV=F, and thus 

uj''''VaH\7bV^ = 2Vc{uj''^yaHybF)VF 

-2Vc(cc;'^Vaff)VbFV^F. (2) 

The first term on the right vanishes because H and F 
commute, whereas the second term vanishes on account 
of the Killing equation satisfied by Z". 

To proceed further we now introduce briefly the ele- 
ments of stochastic differential geometry. The basic pro- 
cess we consider is the Wiener process Wt defined on a 
probability space (il, J^, P), where fi is the sample space, 
J- denotes the filtration, and P is the probability mea- 
sure. We say that Wt is a Wiener process if it is contin- 
uous, Wq = Q,Wt — Ws {Q < s <t) is independent of the 
information Ts up to time s, and Wt — Ws is normally 
distributed with mean zero and variance i — s. A process 
Gt is said to be adapted to the filtration Ft generated by 
Wt if its random value at time t is determined by the 
history of Wt up to that time. If at is ^(-adapted, then 
the stochastic integral Mt = agdWs exists, provided 
at is almost surely square-integrable. If the variance of 
Mt exists, then Mt satisfies the martingale conditions 
E[|Mt|] < oo and E[Mt\Ts] = Ms, where E[-] denotes 
expectation with respect to P. The second condition im- 
plies that, given the history of the Wiener process up to 
time s, the expectation of Mt for t > s is given by its 
value at s. The variance of Mt is determined by the Ito 
isometry E[Mf] = E[J^ a'^ds]. 

A general Ito process is defined by the integral xt = 
Xf) + figds + J*asdWs, where the adapted processes 
fit and at are called the drift and the volatility of xt- A 
convenient way to express this is to write dxt = fitdt + 
atdWt- In the special case fit = fJ-{xt) and at = a{xt), 
where fi{x) and a{x) are prescribed functions, the process 
Xt is said to be a diffusion. 

This analysis can be generalised to the case of a dif- 
fusion Xt taking values on a manifold 9Jl, driven by 
a standard m-dimensional Wiener process {i = 
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1, 2, • ■ • , m). Let Va be a torsion- free connection on 971, 
and suppose ^l°'{x) and cr°(a;) are m -\- \ vector fields 
on 971. Then for the general diffusion on 971 we have 
the stochastic equation dx"" — ii°-dt + afdW^, where dx"" 
is the covariant Ito differential. The quadratic relation 
dx'^dx'' = h''''dt, where h"'' = crfa'", follows from the 
identities dt^ = 0, dtdW^ = 0, and dW^dW^ = 6'^dt. If 
we choose ^i" = 2uj°-^Vi,H, af = 0, and 971 = T, then xt 
reduces to the Schrodinger evolution. 

For any smooth function on 971 we define the 

process (j)t = 4>{xt), and Ito's formula takes the form 
d(j)t — a<\>)dx°- + \{y h<t>)dx°- dx^ , or, more expUcitly, 



dt^alSIaHWl- (3) 



The probability law for xt is characterised by a density 
function p(x, i) which satisfies the Fokker-Planck equa- 
tion dpjdt = -V a^ii" p) + \^ aSI b{h''^ P)- Such a diffusion 
is nondegenerate if h°'^ is of maximal rank. In particu- 
lar, if Qah is a Riemannian metric on 971 and Vq is the 
associated Levi-Civita connection, then if h""^ = cr^g'^'^, 
the process Xt is a Brownian motion with drift on 97t, 
with volatility parameter a. If h"''' is not of maximal 
rank, then the diffusion is degenerate, which is the case 
of relevance to the present consideration. 

Our intention is to generalise the Schrodinger dynam- 
ics to a stochastic process on a nonlinear quantum state 
manifold 971. Specifically, we consider the stochastic re- 
duction model of Hughston ||^ , for which the dynamical 
trajectories are governed by the following stochastic dif- 
ferential equation: 

dx'^ = (2tj°''Vhi/ - icr^V^y^) dt + aV'HdWt. (4) 

When 971 is the state space F of linear quantum mechan- 
ics, then (^) has the following interpretation. The first 
term in the drift generates the unitary part of the evo- 
lution, while the second term creates a tendency for the 
system to evolve to a state of lower energy variance. The 
volatility term is given by the gradient of the Hamilto- 
nian, and generates fluctuations that die down as the 
system approaches an eigenstate. The parameter a con- 
trols the magnitude of the fluctuations. Starting from 
any initial state, the state vector collapses to an energy 
eigenstate, with collapse probability given by the Dirac 
transition probability. Furthermore, if the evolution of 
the density function p{x, t) associated with the process 
is lifted to 7i, we recover the Lindblad form of the 
density matrix dynamics 10|. 

In the case of a general nonlinear quantum phase space 
97t, the stochastic process (|^) can be carried over directly, 
and we obtain the following characterisation of the dy- 
namics of the energy: 

Theorem 1. The Hamiltonian process Ht = H{xt) 
is a martingale, given by the stochastic integral Ht = 
Ho+a V.dWs, where Vt = V"ixt). 

Therefore, under the stochastic dynamics (^), the en- 
ergy of the system is weakly conserved in the sense that 



the expectation of the value of the Hamiltonian at any 
future time is given by the initial value Hq. In particular, 
if reduction occurs, the martingale property ensures that 
the expectation of the terminal value of the energy, which 
is given by the sum of the energy eigenvalues weighted 
by the associated transition probabilities, is Hq. This 
in turn justifies the interpretation of H(x) as the expec- 
tation of the energy in the given state x e 971. The 
proof of Theorem 1 follows from an application of Ito's 
formula (^. The fact that Ht is a martingale does not 
imply that reduction occurs. For a reduction to energy 
eigenstates, we require lim(_,oo Vt = 0. To determine the 
circumstances under which this occurs, we consider the 
dynamics of Vt . 

Lemma 4. The process Vt ~ VaHV^H satisfies 



dVt = a^{VHV''HV''H\7,ybycH)dt 
+<7{\7''H\7aV)dWt. 



(5) 



The proof is as follows. We note that according to Ito's 
formula (H) we have 

dVt = {2uj''''VaVVbH - \a^VaVV^V 

+ \a^V''Hy^HVaybV) dt + (jy'HyaVdWt. (6) 

The first term in the drift vanishes by Lemma 3. The 
remaining two terms in the drift combine to yield (|^), be- 
cause S/'^HV^HVa'^bV = 2V°£rv''7JV'=i?VaVfcVcff + 

^VaVVV. 

For state reduction, we need to show that the drift of V 
is negative. To obtain a suitable criterion for this we pro- 
ceed as follows. If Xa is a Killing field, the cyclic identity 
V[aVfcXc] = impHes that VcVa^b = R^i,fXd, where 
the Riemann tensor is defined by VaVfcAc — VbVa^c = 
—Rabc^d foi' ^'^y vector field Aa- As shown in |jri| we 
thus have: 

Lemma 5. If F{x) is an observable on 971, then 

Next, we define the holomorphic sectional curvature 
ICh of the Kahler manifold 971 with respect to the J- 
invariant plane V^'^Hj''lVH by the formula 



K-H = - 



R, 



apbq 



JPJ\V''HV''HV''HV'^H 



(VaHV-H) 



(7) 



The meaning of K-h is as follows. At each point x G 971 we 
consider the tangent plane spanned by vectors V-ff and 
J^^VH. Then the totality of the geodesic curves tangent 
to this plane at x forms a two dimensional surface in 971, 
and the Gauss curvature of this surface at x is JCh ■ 

Theorem 2. If ICh > 0, then V^ is a supermartingale 
and (^) is a reduction process. 

The proof is by virtue of Lemmas 4 and 5. Writing 
ICt = KHixt), we deduce that 



Vt^Vo- a' 



ICsVfds 



a I VaHVVdWs, (8) 
/o 
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and thus E[Vt|^s] < Vs, the supermartingale condition. 
In particular, if we write Vt = E[T4] then it follows 
from (|)_that dVt/dt = -Ka^Vt^{l + r/f), where Vt^r]t = 
- Vt)^] + K-^E[{JCt ~ K)Vt^], and k = Mm ICh- In- 
tegrating, we obtain Vt = Vb/(1 + Ka^Vo{t + ^t)), where 

— ds. The Hamiltonian sectional curvature is 
positive iff K > 0, in which case > 0. It follows that 
Vt < Vb/(1 + Kcr^Vot), and thus linit^oo Vt — almost 
surely. Therefore, wave function collapse on the nonlin- 
ear state space is guaranteed if the holomorphic sectional 
curvature is positive. The characteristic reduction time- 
scale is r = (kct^Vo)"^, and for t ^ t the uncertainty is 
reduced to a fraction of its initial value. 

Now we are in a position to determine the relation- 
ship between the initial energy uncertainty Vq = V^ (xq) 
and the terminal variance of the energy as a result of 
a reduction. It follows from Theorem 1, together with 
the Ito isometry, that E [{Ht — ffo)^] = o-^Qt, where 
Ha = E[Ht], Qt - E[Qt], and Qt = £ K'^^s. On 
the other hand, if k > then from (0) we obtain 
Vt 1!^ Vq — Ka^Qt- Furthermore, if A = sup^^H and 
A > 0, then (^) implies that Vt >Vo- Xa^Qf Therefore, 
by Theorem 2, if K.^ > 0, we obtain the following bounds 
for the terminal energy variance: 



^> lunE[iHt-Hor]>^. 



(9) 



In particular, if 971 = F, it follows from the relation 

Rabcd = —\ {gacQbd — gbcQad + ^ac^bd — ^bc^-^ad + "^^^ab^cd) 

that Kh = 1 and Vq is the terminal energy dispersion. 

An important issue in the consideration of state re- 
duction processes is whether an energy-based dynamics 
suffices. To address this issue we examine the processes 
induced by for observables other than H. 

Theorem 3. If an observable F commutes with the 
Hamiltonian, then under the stochastic Schrddinger dy- 
namics (^) the process Ft = F{xt) is a martingale. 

The proof follows as a consequence of Ito's lemma with 
an application of Lemmas 2 and 3. Theorem 3 generalises 
a result of [D obtained in the case 9Jt = F. To determine 
whether the system necessarily collapses to an eigenstate 
of F under (0) we require the concept of holomorphic 
bisectional curvature [Q. 

Theorem 4. If the observable F commutes with the 
Hamiltonian, then the stochastic equation for V^ is 



dVt' 



-a^lCFHV^V^dt + aVaHVV^dWt, (10) 



where the holomorphic bisectional curvature of 97t with 
respect to the J-invariant planes determined by F and 
H is defined by 



IC 



FH 



RapbqJ^cJ''d^''F^^H\J''F\J'^H 



(11) 



VaFV''FVbH\/''H 
To prove this result, we use Ito's formula (||) to obtain 
dVt^ = icr^ {y^HV^HVaybV 

- V''V''HVaHVbV^) dt + aVHVaV^dWt. (12) 



Then, by use of Lemmas 2, 3, 5, a calculation shows that 
the two terms in the drift combine to give ([l0|). As a 
consequence we have: 

Theorem 5. If the holomorphic bisectional curvature 
is positive, then for any observable F commuting with 
the Hamiltonian, the associated dispersion V^ is a su- 
permartingale, and (^ is a reduction process for F . 

In the linear theory, this result is intuitively expected, 
because the eigenstates of H also diagonalise any com- 
muting observable F . Indeed, when 9Jl = F we have 
^FH — 5(1 + cos^S), where Q is the angle between 
the vectors V°7? and V°F, from which it follows that 
\ < ICfh < 1, and thus reduction is guaranteed. 
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